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ABSTRACT 

We present here the results of the X-ray analysis of Galactic supernova remnant 
G346.6-0.2 observed with Suzaku. K-shell emission lines of Mg, Si, S, Ca and Fe are 
detected clearly for the first time. Strong emission lines of Si and S imply that X-ray 
emission nature of G346.6-0.2 is ejecta-dominated. The ejecta-dominated emission is 
well fitted with a combined model consisting of thermal plasma in non-equilibrium 
ionization and a non-thermal component, which can be regarded as synchrotron emis- 
sion with a photon index of T ~ 0.6. Absorbing column density of JVh ~ 2.1 X 10 22 
cm -2 is obtained from the best-fitting implying a high-density medium, high electron 
temperature of kT c ^1.2 keV, and ionization timescale of n c t ~ 2.9 x 10 11 cm _3 s 
indicating that this remnant may be far from full ionization equilibrium. The rela- 
tive abundances from the ejecta show that the remnant originates from a Type la 
supernova explosion. 
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1 INTRODUCTION 

A supernova remnant (SNR) consists of expelled material 
called "ejecta" from the explosion and a swept-up interstel- 
lar matter. The X-ray emission results from the interactions 
of ejecta and swept-up matter. From the X-ray emission of 
SNRs we may obtain valuable information about the phys- 
ical properties of the ejecta, swept-up plasma, elemental 
abundances and the history of the explosion. Shock wave 
may be the blast wave associated with the stellar explosion 
and/or the reverse shock wave, which propagates inwards 
from the decelerated blast wave and raises the tempera- 
ture of the stellar ejecta. The young SNRs are bright in 
X-rays and dominated by the emission from ejecta. Thus 
they provide fruitful information about the elements syn- 
thesized by the supernova (SN) explosions. Therefore, the 
observation of the young SNRs is the best method to in- 
vestigate the abundances of the elements synthesized by 
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the SNe (see e.g. SN1006 jKovama et al l Il995l). RCW86 



llBamba. Kovama fe Tomidal bOOcfT Tvcho l|Warren et al.l 
12005ft for the best studied ejecta-dominated SNRs in the 
X-ray band). 

G346.6-0.2 (RA(2000) = 17 h 10 m 19 s , Dec. (2000) = 
—40° 11') is a shell- type SNR located in the Galactic plane. It 
was discovered bv lClark. Caswell fe Green I dl975|) in radio 
band having an angular size of 8 arcmin IjWhiteoak fe GreerJ 
Il99fj ). In X-ray band, on the other hand, G346.6-0.2 was 
first observed by ASC A during its Galactic plane survey 
ijYamauchi et alj|2008h . It is shown that the size of the X-ray 
emission from G346.6-0.2 is less extended than its reported 
radio structure. Five OH(1720 MHz) masers were detected 
toward this SNR and they are all located alo ng the southern 
edge of the remnant (Koralcs kv et al.lll998l ). 

We are studying the X-ray emissions from the ejecta- 
dominated SNRs. For this purpose we have chosen several 
small size Galactic SNRs, one of which is G346.6-0.2. De- 
tailed properties of G346.6-0.2 remained unknown so far be- 
cause of very limited photon statistics in the ASCA (AGPS) 
dat a. Suzaku is the most recent X-ray astronomical satellite 
(see lMitsuda et all (|2007l )) having a large collecting area and 
low background. Therefore, it is the best instrument for ob- 
serving dim and diffuse sources. We propose here to study 
G346.6-0.2 to understand the origin of its thermal and non- 
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thermal X-ray emission which will help us to distinguish its 
SN explosion type. By using the archival data of Suzaku, we 
were able to produce higher quality image and the spectra 
of the remnant, which lead to the results in this study. 

In Section 2, we describe the observation log and the 
data reduction methods. We present the image analysis in 
Section 2.1 and the spectral analysis in Section 2.2. We 
discuss our results and the origin of the thermal and non- 
thermal emission in Section 3. 



2 OBSERVATION AND DATA REDUCTION 

Suzaku satellite has two sets of instruments, one being the 
four X -ray imaging spectrometers (XISs, see iKovama et ah! 
(|2007f )): each of four are at the f ocus of an X-ray tele- 
scope (XRTs, see lSerlemitsos et al. I (120071)) and also a sepa - 
rate hard X-ray detect or (HXD, see lTakahashi et all (|2007h ; 
iKokubun et al. | |2003)). XIS's have two different type of 
CCDs: one being three front illuminated CCDs (namely, FI, 
XISO, XIS2 and XIS3) and other is a back illuminated CCD 
(namely, BI, XIS1). Each CCD covers an area of 17.8x17.8 
arcmin 2 . 

The data used in this analysis is taken by the XIS on- 
board Suzaku. G346.6-0.2 was observed on 2009 October 07 
for 56.7 ks (Obs ID:504096010). During the observation, the 
XIS were operated in the normal full-frame clocking mode 
with the editing mode of 3 x 3 and 5x5, which are low and 
medium data rates, high and super-high data rates, respec- 
tively. 

For the data reduction and spectral analysis, we used 
the HEADAS soft ware package o f version 6.5 and xspec ver- 
sion 11.3.2 (see[A rnau 1 (|l996l )). The response matrix files 
(RMF) and ancillary response files (ARF) were made us- 
ing xisrmfgen and xissimarfgen version 2006-10-17 (see 
llshisaki et all (120071 )). 



2.1 Image Analysis 

Figure 1 shows the XISO image of G346.6-0.2 in the 0.3 - 
10 keV energy band. The solid circle (radius ~5 arcmin) 
from the centre) shows the region where the spectra is ex- 
tracted. The background region is shown by dotted circle 
(RA(2000) = 17 h 10 m 54 s , Dec. (2000) = -40°06'24", radius 
~1.8 arcmin). We excluded the upper right corner of the de- 
tectors containing the onboard masked out 55 Fe calibration 
sources. For comp arison radio continuum im age of G346.6- 
0.2 at 843 MHz bv lWhiteoak fe Greenl (|l996l ) is overlaid on 
the figure, where four maser sources a r e also pointed out as 
has been reported by iKoraleskv et all (|l998h . 

The bright image in XISO is clearly seen and we have 
obtained the Si-S map in 1.5 — 2.5 keV energy band (Figure 
2) to see if this emission is from interstellar medium (ISM) 
or from the ejecta. In Fig. 2, Si-S emission is distributed 
throughout the remnant and it is bright, which may well 
indicate that the X-ray emission coming from the remnant 
is ejecta dominated rather than originating from the ISM. 



2.2 Spectral Analysis 

The extracted XIS spectra from a circular region of G346.6- 
0.2 (see Fig. 1) are analyzed in the full energy band of 0.3— 10 



keV. Figure 3 shows representative spectra of XISO, XIS1 
and XIS3 simultaneously, extracted from the region (shown 
in Fig. 1 by the solid circle) with its corresponding best- 
fitting model and residuals. 

The strong Si and S lines from the spectra and Si-S 
image (see Fig. 2) in 1.5 — 2.5 keV energy band suggest that 
the emission is ejecta dominated. Therefore, we have applied 
one-component models in our spe ctral fits. We first fitted 
the da ta with an absorbed (wabs; iMorrison fe McCammor] 
(1983)) VNEI, which is a model for a non-equilibrium ion- 
ization (NEI) collisional plasma with variable abundances 
|Borkowski. Lverlv fe Reynolds! [200 ll ). The model did not 
fit well (reduced x 2 of 1.12 (637/567 d.o.f.)). So, we added 
a power-law model (non-thermal synchrotron emission from 
relativistic electrons). The parameters of the absorbing 
column density (Ah); electron temperature (fcT c ); the 
ionization parameter (n e t) and the metal abundances C, 
N, O, Ne, Mg, Si, S, Ar, Ca, and Fe are set free in this fit 
resulting in an unacceptable error value. Then, we set Mg, 
S, Si, Ca, and Fe free, since their emission lines are clearly 
seen in the spectra, while the rest are set to solar values 
jAnders fe Grevessel Il989l ) resulting in a reasonably good 
reduced \ 2 value of 1.07 (606/565 d.o.f.) with reasonable 
error limits. Similar steps were also applied to VMEKAL 
model, which describes an emission spectrum from hot 
ionized gas in collisio nal ionization equilibrium with vari- 
able abundances dMewe. Gronenschild fe van den Oordl 



Mewe. Lemen fe van den Oordl 1 19861 ; 

Liedahl. Osterheld fe Goldstein I Il995l ). and VPSHOCK 



19851; 



model, which is suitable for modelling plane-parallel shocks 
in young SN Rs, where plasma has not reached the io nization 
equilibrium (|Borkowski. Lverlv fe Reynolds! l200ll ). These 
two models gave us reduced \ 2 value of 1.22 (691/568 d.o.f.) 
and 1.11 (629/567 d.o.f.), respectively. Then we added a 
power-law component to both models, which yielded almost 
similar results (of reduced xVd.o.f. value of 607/566 and 
605/565, respectively, of the best-fitting parameters and 
errors values) as VNEI model as given in Table 1. 

These three models basically represent the emission 
from optically thin thermal plasma with slight differences. 
For example, VNEI is a NEI model similar to VPSHOCK, 
except with a single ionization timescale. 

From an absorbed VNEI and power-law model, the ab- 
sorbing column density, the electron temperature, the ion- 
ization parameter, the metal abundances of Mg, Si, S, Ca, 
Fe, volume emission measure (VEM) and flux (Fx) of the 
VNEI component and the photon index (r) and norm of 
power- law component are obtained and given in Table 1, 
with the corresponding error values with 90 per cent confi- 
dence levels. Total flux and reduced x 2 /d.o.f. value of VNEI 
and power-law components are also given in this table. The 
results of VMEKAL and power-law, VPSHOCK and power- 
law models are also presented in Table 1 for comparison. 

To determined the line centre energy of the K- 
shell (Ka) lines, we fitted the spectra using a simple 
bremsstrahlung continuum plus five Gaussian lines with an 
absorbing column. The obtained best-fitting central ener- 
gies of Ka emission lines are He-like Mg (1.38±0.06 keV), 
Si (1.86±0.02 keV), S (2.52±0.06 keV), Ca (3.25±0.55 keV) 
and Fe (6.77±0.32 keV). 

Figure 4 shows the best-fitting met al abundances (nor- 
malized to Si) relative to solar values l|Anders fe Grevessel 
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1 19891 ) with the predicted n ucleosynthesis W7 Type la SN 
model (|Nomoto et al.l ll997T ) , which is widely used as the 
standard model. 



3 DISCUSSION AND CONCLUSIONS 

In this work, we provide for the first time a high-quality 
image, spectra, and detailed description of the X-ray emis- 
sion of G346.6-0.2 using a set of public archival Suzaku 
data. The radio emission of th is remnant extends to 8 ar- 
cmin (Whitcoa k fc Greenlll996T ). while X-ray image shows a 
smaller extension (see Fig. 1), which is consistent wit h that 
of the ASCA results shown in lYamauchi et al.l (|2008T ). 

We fitted the spectra with "thermal" component in 
non-equilibrium ionization and an additional "non-thermal" 
power-law component with reasonably good chi-square val- 
ues. 

Thermal emission. The thermal part of the spectra can 
be represented by non-equilibrium ionization (VNEI) model. 
The thermal X-ray emission of young SNRs, whose emission 
is still dominated by the ejecta, is predominantly from the 
ejecta heated by the reverse shock. Since the ejecta is metal 
abundant, its X-ray spectra usually shows strong emission 
lines of heavy elements such as O, Ne, Mg, Si, S, Ar, Ca and 
Fe. Ka emission lines of Mg, Si, S, Ca and Fe for G346.6-0.2 
are detected clearly for the first time in this work. 

The interstellar absorbing column density JVh is ob- 
tained to be (2.1±0.2)xl0 22 cm -2 , which is consistent with 
the galactic value of (1.5± 0.2)xl0 22 cm -2 in that direction 
jDickev fc Lockmanlll990l ). ASCA r esults of G346.6-0.2 re - 
ported a value of N H > 10 22 cm" 2 (lYamauchi et al.ll2008h . 
indicating that X-ray spectra of the remnant is heavily ab- 
sorbed by interstellar matter, perhaps due to its location 
being close to the Galactic plane. In other words, the elec- 
tron density (n e ) of medium of the remnant is high enough. 
This may also be supported by the fact that the south region 
of the remnant is in interaction with f our m olecular clouds 
as has been noted by iKoraleskv et all l|l998t) . 

Being close to the Galactic plane and being in a highly 
dense region, one can expect that some portions of the rem- 
nant can be highly asymmetrical. The asymmetry in the 
south region of the remnant may well be due to the inter- 
action with the maser sources as has been noted in Fig. 1. 
One can argue that the surface brightness of the remnant is 
far from the equal distribution and spherical symmetry due 
to the inhomogeneity of the medium. 

The thermal emission further provides us information 
about the electron temperature and age of the plasma. Our 
best-fitting temperature is obtained to be ~1.2 keV, which is 
typical for shell-like SNRs. This value is also consistent with 
the result of ~1.6 keV obtained from ASCA observations 
lYamauchi et alJl200Sh . 

The age of G346.6-0.2 has not been predicted so far. Es- 
pecially, young SNRs whose ages are a few hundred or about 
a thousand years are not only bright but also dominated by 
the emission from ejecta. Being an ejecta-dominant rem- 
nant and bright in X-rays G346.6-0.2 can be a young SNR. 
Plasma in most of young SNRs are in the NEI condition. The 
ionization age (or ionization parameter) is denned as r—n c t, 
which is often used as a key diagnostic of the NEI stat e, r 
is typically required to be ^10 cm" " 3 s l|Masailll984l ) for 



full ionization equilibrium. In our observations we obtained 
n e t to be (2.92±0.01) x 10 11 cm _3 s, which may also indicate 
that this SNR is far from full ionization equilibrium. So, 
the plasma has not yet had time to reach ionization equi- 
librium, and it is still being ionized. We calculated n e to be 
~0.82 cm -3 from VEM=n e nnV where n B = 1.2riH and V is 
the X-ray-emitting volume of the remnant (estimated to be 
~2 x 10 54 cm 3 ) by adopting an average distances of d ~ 8.3 
kpc (average of a near value of 5 .5 kp c and a far value of 
11 kpc given by IKoraleskv et all (|l998h ) and assuming the 
emitting volume to be spherical shell. From the equation 
t=r/n e , the age of G346.6-0.2 calculated to be ~1.1 x 10 4 
yr. Thus, this remnant is likely to be the oldest known ejecta- 
dominated shell-like SNR ( e.g. SN1006, RCW86, G3 37.2-0.7 
an d G309. 2-0.6, ~1000 yr ilYamaguchi et al.ll2008ah . ~1800 
yr (lYamaguchi et alj|2008bh. 2000 -4500 vr and 700-4000 yr 
(Rak owski. Hughes fc Sl anc 2001), respectively.) 

Non-thermal e mission. Recently, in some SNRs 
(e.g. first S N1006 llKovama et ail 119951) and then RX 
J1713.7.3946 dKovama et all Il997l; Isiane et aljfl999h . Cas 



A dVink et~alT 2000l ; iHughes et al.N2000D . and Tvcho's SNR 
l|Hwang et al.l [2002)) non-thermal emission has been de- 
tected. There are two possible mechanisms for non-thermal 
emission: one is synchrotron emission and the another is non- 
thermal bremsstrahlung. The non-thermal emission com- 
ing from G346.6-0.2 is most likely a synchrotron emission 
(power-law component with photon index of F ~ 0.6). The 
real source of synchrotron emission in SNRs is believed to 
be the acceleration of the ultrarelativistic electrons in the 
magnetic fields with the help of the shock waves. A power- 
ful X-ray point source inside the remnant or in the vicin- 
ity of the remnant can also contribute to the non-thermal 
emission. The photon index we obtained (F ~ 0.6) is lower 
(harder) than those of typical shell-like SNRs. This forces us 
to presume that a pulsar wind nebula (PWN) contribution 
is also possible, although a PWN has not been discovered so 
far in or in the vicinity of G346.6-0.2. 

The flux value of G346.6-0.2 is obtained to be F x ~ 4.4 x 
10 -11 ergs _1 cm -2 for 0.3 — 10 keV range, which is in a good 
agrement with that of found with ASCA (|Yamauchi et al.l 
I2008T I. Using this flux value for the source and also taking 
its X-ray angular size of ~5 arcmin, its surface brightness in 
the 0.3 - 10 keV energy range is found to be E ~ 1.8 x 10" 12 
erg s _1 cm _2 arcmin -2 , which is higher than many shell-like 
Galactic SNRs. Surface brightness mainly depends on the 
explosion energy and the density of the medium, but a PWN 
may also contribute to the surface brightness of the SNR. In 
this case, high medium density and small size of the remnant 
may be the main reason of high X-ray surface brightness of 
G346.6-0.2. 

Relative abundances in the ejecta. Theoreti- 
cal models predict that in core collapse SN explo- 
sion low-Z elements like O, Ne, Mg are produced 
(I Thieleman n. Nomoto & Hashimotol I1996T) . while in Type 
la SNe (jNomoto et alj ll99Tt > high-Z elements like Ar, Ca 
and Fe are mostly produced. Furthermore, Fe production 
in Type la SNe is far larger than that of core-collapse SNe 
|Nomoto et al.lll997l ; llwamoto et al.lll999T ). The detection of 
Fe emission lines from the remnants is important to identify 
the type of the remnant. Suzaku has detected Fe Ka lines 
in especially Type la SNRs (e.g. Tycho, SN1006, RCW86). 

When we compare our best-fitting abundances of Mg, 
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S, Ca and Fe relative to Si with the predicted theoret- 
ical values, Mg relative to Si is higher, while S and Ca 
relative to Si are lower. The reason for the lower abun- 
dances of S and Ca could be because the elements are 
concentrated at the inner layers of the remnant and hence 
are not heated enough by the reverse shock yet. This in- 
formation implies that this SNR is at its early evolution 
phase. Our best-fitting abundance Fe relative to Si is con- 
sistent with the expected value (see Fig. 4) within the con- 
fidence range. However, considering the large error bars 
due to low statistics, we could also comment a similar 
statement for Fe abundance value and an unattained re- 
verse sh ock, as it is reported in many T ype la SNRs like 
SN10 06 <|Yamaguchi et al.ll2008al). Tvcho l|Tamagawa et al] 
I2009D and G337.2-0.7 l|Rakowski. Hughes fc Slaneil200ll ) "the 
Fe abundance is being less than expected value produced in 
a Type la SNe. Considering all these cases we can predict 
that the remnant may be originated from a Type la SN ex- 
plosion. 
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Table 1. The best-fitting parameters of the G346.6-0.2 spectra in 0.3 — 10 keV energy band for VNEI, 
VMEKAL and VPSHOCK models. 



Parameters 


VNEI+power-law 


VMEKAL±power-law 


VPSHOCK±power-law 


7V H (xl0 22 cm" 2 ) 


2.1 ±0.2 


2.0 ±0.2 


2.1 ±0.2 


fcT c (keV) 


1.22±0.04 


0.97±0.05 


1.3±0.1 


n c t(xl0 1:L cm _3 s) 


2.92±0.01 




(6.92±0.01^, 


Mg 


0.51±0.12 


0.65±0.15 


0.76±0.16 


Si 


0.47±0.07 


0.44±0.07 


0.56±0.09 


S 


0.7±0.1 


0.7±0.1 


0.8±0.1 


Ca 


2.3±0.6 


3.4±0.8 


2.3±0.7 


Fe 


0.6±0.3 


0.4±0.2 


0.6±0.3 


VEM^ 


11.6±1.1 


15.1±1.9 


10.4±1.5 


FluxJ 


3.1±0.5 


2.9±0.6 


3.2±0.4 


Photon Index 


0.6 ±0.3 


0.5 ±0.3 


0.6 ±0.3 


norm(xl0 — 4 photons cm~ 2 s -1 ) 


4.46 ± 0.04 


4.55±0.04 


4.55±0.05 


Flux£ 


4.4±0.2 


4.5±0.2 


4.4±0.2 


X 2 /d.~f. 


606/565=1.07 


607/566=1.07 


605/565=1.07 



a The VPSHOCK fitting gives an upper limit on the ionization timescale. 

Volume emission measure (VEM= J n c njidV in the unit of 10 cm , where n c and rtn are number 
densities of electrons and protons, respectively and V is the X-ray-emitting volume. 
c Flux in the 0.3 — 10 keV band in the unit of 10 -11 erg s _1 cm -2 . 
^ Total flux in the 0.3 — 10 keV band in the unit of 10 -11 erg s _1 cm -2 . 




counts/pixel 2 4 6 8 10 12 

Figure 1. XIS0 image of G346. 6-0.2 in the 0.3 - 10 keV full 
energy band. The coordinates are referred to epoch J2000. The 
calibration source located at the upper right corner, as seen by 
red-circle, is masked out. The colour-coding is for the X-ray pho- 
tons in counts/pixel. The overlaid iso-intensity contours are from 
radio observation at 843 MHz. The contour levels are -7, 6, 25, 50, 
95, 135, 165, 190, 210, 230, 255 and 270 mjy/bcam, selected the 
same as the radio works for visual comparison. The regions used 
to extract spectrum and to determine the background parameters 
are indicated by blue solid and dotte d circles, respectively. The 
positions of the four masers defined by Koralesky et al. ( 1998]) are 
marked as crosses. 
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Figure 2. G346.6-0.2 XIS0 image in 1.5 - 2.5 keV (Si-S) energy 
band. The image is smoothed for 10cr Gaussian in order highlight 
the Si-S distribution. Overlaid contours are spaced linearly in 
intensity (0.01, 0.02, 0.03, 0.05, 0.09 and 0.16 counts/pixel). 
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Figure 3. Background-subtracted XIS (XISO:black, XISl:red, XIS3:green) spectra of G346.6-0.2 in full energy band (0.3-10 keV) fitted 
with an absorbed VNEI and power-law model. The lower panel shows the residuals from the best-fitting model. 
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Figure 4. Metal abundances of Mg, S, Ca and Fe relative to solar 
values and normalized to Si. W7 model for a Type la is shown by 
diamonds while our data are shown by triangles. 



